import numpy as np

# 路面结构强度指数（PSSI）
# SSR —— 路面结构强度系数（Pavement Structure Strength Ratio），为路面容许弯沉
# 与路面实测代表弯沉之比；
# lR—— 路面容许弯沉（mm），路面容许弯沉计算方法见附录 C；
# l0—— 路面实测代表弯沉（mm）；
# a0—— 模型参数，采用 15.71；
# a1—— 模型参数，采用-5.19。

L0 = 6
Lr = np.random.randint(1, 10)


def PSSI():
    a0 = 15.71
    a1 = -5.19
    ssr = L0 / Lr
    return 100. / (1 + a0 * np.exp(a1 * ssr))


if __name__ == '__main__':
    pssi = PSSI()
    print(pssi)